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We investigate witii a microscopic model anharmonic K-cation oscillation observed by neu- 
tron experiments in /3-pyroclilore superconductor KOS2O6, which also shows a mysterious 
first-order structural transition at Tp— 7.5 K. We have identified a set of microscopic model 
parameters that successfully reproduce the observed temperature dependence and the super- 
conducting transition temperature. Considering changes in the parameters at Tp, we can explain 
puzzling experimental results about electron-phonon coupling and neutron data. Our analysis 
demonstrates that the first-order transition is multipolar transition driven by the octupolar 
component of K-cation oscillations. The octupole moment does not change the symmetry and 
is characteristic to noncentrosymmetric K-cation potential. 
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The /3-pyrochlore oxides A0s20e{A=K, Rb, or Cs) are 
one of the material families that have unique cage-like 
structure, and ^-cation inclusions exhibit large anhar- 
monic local oscillations. '^^^ Such anharmonic and large- 
amplitude oscillations strongly interact with conduction 
electrons, which leads to strong coupling superconduc- 
tivity, ^'^"^^ anomalous temperature dependence of elec- 
tric resistivity'*' and nuclear relaxation time.^' To ex- 
plain such interesting properties is a big challenge for the 
theory of electron-phonon systems and strong coupling 
theory of superconductivity.^^ 

As the ^-cation size decreases, amplitude of the A- 
cation oscillation increases, and thus potassium oscilla- 
tions are the most anharmonic. Elastic neutron scat- 
tering data show that the temperature dependence curve 
of oscillation amplitude is concave (convex upwards) in 
KOS2O6, reflecting strong anharmonicity.*^' A recent 
inelastic neutron scattering experiment revealed that the 
phonon peaks shift to lower energy with decreasing tem- 
perature and the peak positions are at around 3-7 meV at 
the lowest temperature.*^ This softening is the strongest 
in K-, and the weakest in Cs-compound. Superconduct- 
ing transition takes places in all the three, and the tran- 
sition temperature is Tc —9.6, 6.3 and 3.3 K for A=K, 
Rb and Cs, respectively.^'^ The symmetry of the super- 
conducting gap function is confirmed to be s-wave.^'^ 

KOS2O6, the most anharmonic one, not only has the 
highest Tc of superconductivity, but also exhibits an- 
other singularity, a first-order structural transition at 
Tp — 7.5 K.^ It is interesting that no sign of symmetry 
breaking has been observed.*' ^°'^* Below Tp, the tem- 
perature dependence of electric resistivity changes to T^ 
from r''(7 ^ 0.5) at higher temperatures,* and the abso- 
lute value decreases by about 25 % at Tp. Both of these 
results show that electron-phonon scatterings are sup- 
pressed below Tp,^ which is also supported by the reduc- 
tion in the specific heat jump at Tc in magnetic fields 
H when Tc{H) < Tp{H) ~ Tp{0).'^ However, recent neu- 



tron scattering data^° show the oscillation amplitude in- 
creases below Tp. Its naive interpretation is an increase 
in the electron-phonon coupling. 

In our previous work,^^ we have developed a general 
theory of anharmonic ion oscillations in tetrahedral sym- 
metry. A strong coupling theory of superconductivity 
mediated by these ion oscillations has been also devel- 
oped. Applying them to /^-pyrochlore compounds, we es- 
timated their T^. We also started to study the contradic- 
tion between the change of the electron phonon coupling 
and the neutron data at Tp in KOS2O6 but that was in 
a qualitative level. 

In this Letter, we will resolve by microscopic calcu- 
lations the contradiction between the data of resistiv- 
ity and specific heat'* and the neutron scattering^" in 
KOS2O6 at Tp. Our calculation demonstrates that the 
first-order transition is a multipolar isomorphic phase 
transition and the phonon anharmonicity plays a cru- 
cial role to explain the experimental data. We will apply 
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Fig. 1. (color online) Potential profile, and probability density 
of the ground state wavefunction |\I'ga(r)[^ for K (solid), Rb 
(dashed), and Cs (dotted line), respectively, along [111] direc- 
tion. Energy eigenvalues are indicated by short bars and the 
numbers represent the degeneracy. 
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the theory developed in ref. 22 to AOS2O6 and carry out 
more elaborate numerical calculations in order to under- 
stand the higher-temperature and higher-energy prop- 
erties of the anharmonic ^-cation oscillations. We will 
show that we can reproduce the temperature dependence 
of the softening of phonon energy and the amplitude of 
A-cation oscillations in a wide range of temperatures in 
nice agreement with the experimental data. 

We start with a short review of our model for A-cation 
oscillation. Since the inelastic neutron experiment shows 
nearly momentum-independent modes of their dynam- 
ics/^ we employ a local model. The A-site has the tetra- 
hedral Td symmetry and this implies that the ion po- 
tential generally has, in addition to spherical and cu- 
bic fourth order terms, a third-order anharmonic term, 
which breaks space inversion symmetry. The ionic Hamil- 
tonian is thus given by 

V72 71, ,r 2 

H^on = -^ + ^\r\'+bxyz + c,\r\* + C2?%l) 

where h is set to be unity, and r — (x, y, z) is the ion dis- 
placement from the equilibrium position, = x'^+y'^+z'^ 
and M is the ion mass. We ignore the higher order terms 
of 0(|r|^) which is irrelevant for our discussions because 
of positive forth order terms. We diagonalize Hamilto- 
nian (1) numerically in the restricted Hilbert space of di- 
mension 62196, which is about three times larger Hilbert 
space used in the previous study^^ and sufficient to dis- 
cuss the physical properties below room temperature, by 
extrapolating the obtained data to the infinite limit if 
necessary. Details of calculations are explained in ref. 22. 

An important point is the variation in potential pa- 
rameters among different compounds. The Madelung 
part of the potential energy is essentially the same among 
the three compounds, since the band structure calcula- 
tions^^ show essentially the same electric states for them. 
Thus, the main difference in the potential parameters 
originates from the the relaxation of local charge density 
as discussed in ref. 23. This point can be well captured 
by setting the smaller uj for the smaller ion. We use the 
same parameters as those in ref. 22: uj = 26.4 K, 54.6 K 
and 74.8 K for K, Rb, and Cs, respectively, with keeping 
the same values for other parameters, b — 9324 K/A'^ 
and ci = 4c2 — 3332 K/A^. The validity of our choice 
will be checked later by comparing the phonon energy 
and oscillation amplitudes with experimental data. 

Figure 1 shows the potential profile with these pa- 
rameters along [111] direction for the three compounds. 
Eigenenergies of Hamiltonian (1) are indicated by bars 
and the probability density of the ground state |\E'gs(r)|^ 
is also shown. The potential of K-cation is much shal- 
lower than those of Rb and Cs, and correspondingly, its 
|^'gs(r)p has a broader tail than the others. Kunes et at, 
proposed a potential with uj^ < for KOS2O6 based on 
the first principle calculation,^"^ but with their choice, the 
K-cation oscillates about 1 A, which is too large and in- 
consistent with the neutron experiment data.^^' There- 
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Fig. 2. (color online) Temperature dependence of the variance of 
ion oscillation {x^). Three solid lines represent calculated {x'^) 
for the three compounds. Symbols and dotted line show the ex- 
perimental data taken from refs. 16, 17 and 18. 

fore, we use positive for all the three compounds. As 
we will show later, our potential parameters reproduce 
the correct temperature dependence of the variance of 
ion oscillation (x^), and also the excitation energy agree- 
ing with the neutron and x-ray scattering data.-'^^^-'^^ 

We first calculated the variance of ion oscillation {x'^) 
as a function of temperature from the calculated eigen- 
functions of (1) and the results are shown in Fig. 2. The 
experimental values determined by neutron^^' and x- 
ray^^' scattering are also shown for comparison. As 
is clearly seen, our results are quantitatively consistent 
with the experimental results. For KOS2O6, we mainly 
concentrate on fitting the recent experimental data of ref. 
18 and the low-temperature limit is consistent with the 
recent data ^ 0.014 A"^ just above Tp.^° As for Rb- and 
Cs-compounds, the calculated (a;^) is slightly larger than 
the data in ref. 17. However, concerning the discrepancy 
between the data in refs. 17 and 18 for K-compound, 
our results for Rb- and Cs-compounds well capture the 
essential aspect of the experimental data. We note that 
the values in the low-T limit are important to discuss the 
superconducting transition temperature, since these val- 
ues directly infiuence the dimensionless electron-phonon 
coupling constant A. 

Figure 3 shows the phonon spectrum, i.e., the imag- 
inary part of the phonon Green's function D{i>) = 



^iut—rjt 



x(0)]) at six different tempera- 



tures for A—\^, Rb and Cs compounds, r/ is a phenomeno- 
logical relaxation rate set as ry = 3.5 K. There are several 
peaks in the spectrum of KOS2O6 for intermediate tem- 
peratures in Fig. 3(a) owing to the anharmonic terms 
and the smaller u compared to other members. The re- 
sults are consistent with the inelastic neutron scattering 
spectra, particularly in the following three aspects, (i) 
The peak positions at the lowest temperature are 38.5, 
52.9 and 73.1 K, which correspond to the energy of the 
first excited states A, consistent with the neutron data 
at 1.5 K:i9 n"^^ = 38, 56, and 66 K for K, Rb and Cs, 
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Fig. 3. (color online) Phonon spectrum vs energy at six different 
temperatures for (a) K, (b) Rb, and (c) Cs compounds. Spec- 
trum is broaden owing to the many Lorentzians whose width is 
given by 77 = 3.5 K and this is the direct consequence of the an- 
harmonicity. Inset: Temperature dependence of the highest peak 
position for K (circle), Rb(square) and Cs(diamond) taken from 
ref. 19. Lines represent the present calculations. 

respectively, (ii) The softening of the peak energy with 
decreasing T is the strongest for K, intermediate for Rb 
and the weakest for Cs. The temperature dependence is 
also qualitatively reproduced by our calculation as shown 
in the inset of Fig. 3. Note that the relevant mode splits 
into two branches owing to two cages in the unit cell, and 
our data correspond to their average energy, (iii) The 
spectrum of K at higher temperatures is much broader 
than that of Rb and Cs. 

As discussed above, our model can explain the neu- 
tron experiments and well describes potential of A-cation 
oscillations. These results confirm our previous work,^^ 
in which the superconducting transition temperature 
was calculated with the same parameters used in this 
Letter and the result is Tc = 10.5, 5.7 and 3.4 K, for K, 
Rb and Cs, respectively. The calculated Tc agrees with 
the experimental data within ±1 K, even without opti- 
mizing parameters for each compound. The variation in 
Tc among the three members is mainly attributed to the 
difference in the anharmonicity of A-cation oscillations. 
Indeed, 60 % of Tc is owing to the K-cation oscillation in 
KOsaOg. 

As mentioned in the introduction, KOS2O6 exhibits 
a first-order structural transition at Tp = 7.5 K, which 
does not break any symmetry.''' ^'^"'^'^ In ref. 24, we have 
proposed a scenario of isomorphic structural transition to 
explain it based on a simple toy model. The point is that 
in the tetrahedral symmetry, not only (r^) = 3(a;^) = 
3u2 but also the third moment \{xyz)\ = U3 are nonzero 
and their changes at Tp do not break the point group 
symmetry. 

Here, we will examine this scenario by means of the 
present realistic model. Analyzing the changes in the os- 
cillation profile as a function of variations in the poten- 
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Fig. 4. (color online) Contour plot of (a) STc in uj'-b' plane 
with fixing Cj=5c2=ci, (b) SX, (c) i5A, (d) SU2/U2, and 
(e) SU3/U3. In (a), all the five regions are indicated; A: 
(sgn((5A), sgn(5A), sgn(5«2), sgn{5us)) = {+,-,+,+), 
B: {+,-,+,-), C: {+,+,+,-), D: {-,+,+,-) and E: 
(—,-!-,—,—). The solid lines with symbols represent the bound- 
aries and the dashed lines represent the contours. Region D is 
indicated by shade. 

tial parameters at Tp, we will identify the transition as 
"multipolar'^ isomorphic one, because the change in the 
octupole U3 is much larger than in the isotropic scalar 
U2. This point is important to explain the experimental 
results of recent neutron scattering^" and specific heat 
jump at Tc,^ which will be explained below. Our previ- 
ous analysis based on the simplified toy model^"* fails to 
describe this point, since the model cannot distinguish 
U2 and W3 owing to the lack of the degrees of freedom. 
This kind of phase transition driven by a scalar order 
parameter is similar to that proposed for the f-electron 
system PrRu4Pi2-^^ 

Let us discuss the nature of the transition at Tp and 
three experimental results are important. First, the mag- 
netic field dependence of the specific heat jump at Tc"* 
indicates that the electron-phonon coupling A decreases 
below Tp. Secondly, a recent neutron experiment shows 
that U2 jumps up at Tp and the value is ~ 0.02 at 
low temperature 1.5 K.^" Thirdly, the upper critical field 
Hc2{T) shows that the first-order transition does not af- 
fect Tc{H — > 0).'''^® The first two experimental data can- 
not be explained on the basis of harmonic oscillations. In 
the harmonic case, A oc and thus, the increase in U2 
leads to enhanced A. These puzzling experimental results 
are naturally explained on the basis of the multipolar iso- 
morphic transition. 

Since the transition is first-order, the potential pa- 
rameters change discontinuously, which originates from, 
for example, the changes in volume, oxygen positions, 
and the inter-site ion interactions. The high-T poten- 
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tial parameters w==26.4 K, 6=9324 K/A^, ci=4c2= 3332 
K/A'' change to low-T ones, say, w', b', c[ and C2, respec- 
tively. Changes in b and ui are particularly important and 
we investigate the effects of their change, while we set 
c[ =5c2= ci for simplicity. Thus, with varying lj' and b', 
we repeat the same procedure as before^^ and calculate 
the deviations of dimensionless electron-phonon coupling 
SX=X'—X, and gap SA, similarly, Su2, Su^, and ST^.. Here, 
the quantities denoted with prime symbol are the low-T 
values calculated for the new potential parameters at 
for each of the parameter set. Tc is low enough compared 
with A and thus these quantities are regarded as those 
for the low-temperature limit. The results are shown in 
Fig. 4. 

First, it is noted that there exist five regions, A-E, 
distinguished by the sign of i5A, 6 A, 6u2 and SU3. In the 
harmonic case, only the variations characterized by the 
region A or E are possible, but they are not consistent 
with the experimental results, (5A < and 6112 > 0. The 
other three regions appear as a consequence of the third- 
order phonon anharmonicity b. Among them, it is the 
region D that satisfies the experimental constraints, and 
the low-T parameters b' and w' should be located inside 
this region. The change in Tc is less than 2 % and this 
small change is consistent with the last one of the three 
experimental points mentioned before. 

The puzzling behavior, < and 5u2 > 0, is a con- 
sequence of two competing effects. The reduction in uj 
enhances A, U2, and U3, whereas the reduction in |6| sup- 
presses them. The change in A is dominated by the effect 
of reduced while the change in U2 is rather due to 
reduced co. Our calculation shows that the experimental 
results predict the constraint on the potential changes at 
Tp as 2Su}/uj < Sb/b < 1.56io/oj < 0, as far as the changes 
are not so large. 

One can also explain the reduction in the electrical re- 
sistivity observed at Tp,'* since the changes SA > and 
(5A < mean the suppression of the electron scatterings 
due to the K-cation oscillations. Finally, the nature of the 
first-order transition is identified by examining various 
moments of ion oscillations. As for the changes in U2 and 
Us, the isotropic part (r^) = 3u2 of the oscillation am- 
plitude is enhanced while anisotropic one \{xyz)\ = (/;; is 
suppressed as shown in Figs. 4(d) and (e). Furthermore, 
the relative change in us is much larger than that in U2, 
and this is the reason to call this a multipolar isomorphic 
transition. 

The important point is that our calculation already 
reproduces changes consistent with the experimental re- 
sults at Tp. In order to carry out the further quantitative 
comparison between the theory and the experimental re- 
sults, one needs more detailed information about the po- 
tential changes. It is also useful to observe the anisotropic 
part of the Debye- Waller factor and compare with theo- 
retical calculation. 

In summary, we have proposed that the first-order 
transition observed in KOS2O6 is a multipolar isomor- 



phic transition and demonstrated that the third order 
fluctuation (xyz) is the primary order parameter. Our 
theory naturally explains both the reduction in electron- 
phonon coupling constant and the enhanced oscillation 
amplitude at Tp. We hope detailed neutron scattering 
experiments detect its change at Tp. 
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